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Abstract 

We calculate nonperturbative 0(A 2 Q CD /m 2 ) corrections to the dilepton invariant mass spec- 
trum and the forward-backward charge asymmetry in B — > X s e + e~ decay using a heavy quark 
expansion approach. The method has recently been used to estimate long-distance effects in 
B — ► Xsj. We generalize this analysis to the case of nonvanishing photon invariant mass, 
q 2 7^ 0, relevant for the rare decay mode B — > X s e + e~ . In the phenomenologically interesting 
q 2 region away from the cc resonances, the heavy quark expansion approach should provide 
a reasonable description of possible nonperturbative corrections. In particular this picture is 
preferable to the model-dependent approach relying on the tails of Breit-Wigner resonances, 
which has been employed so far in the literature to account for these effects. We find that 
the 0(AQ CD /m 2 ) corrections to the dilepton invariant mass spectrum and to the forward- 
backward asymmetry in B — > X s e + e~ amount to several percent at most for q 2 /m 2 ^0.3 and 
q 2 /m 2 ^ 0.6. The 0(Aq CD /m 2 .) correction to the B — > X s ~f decay rate is also computed and 
found to be +3%, which agrees in magnitude with previous calculations. Finally, we comment 
on long-distance effects in B — > X s vv, which in this case are extremely suppressed due to the 
absence of virtual photon contributions. 
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97-2410, the Monell Foundation and the Seoam Foundation Fellowships. 



1 Introduction 



Inclusive rare decays of B mesons, such as B — > X s 7, B — > X s e + e~ or 5 — > X s vi>, provide 
important opportunities to test the Standard Model of flavor physics. These processes are 
particularly suited for this purpose, since they occur only at the loop level and are dominated 
by contributions that are reliably calculable in perturbation theory. Heavy quark (l/m&) 
expansion and renormalization group improved perturbative QCD form a solid theoretical 
framework to describe these decays. The rates of the B meson decays are essentially deter- 
mined by those of free b quarks, which have been computed at next-to-leading order (NLO) 
in QCD for B -> X sl g §, B -> X s e + e~ [g, § and 5 -> X s z/z/ §. The first subleading 
C?(l/ra&) power corrections in the heavy quark expansion have also been studied and are in 
general well under control || 0, ||, [|. However, to properly assess the prospects for precise 
tests of the Standard Model and its extensions (see for instance [[HJ), further possible sources 
of theoretical uncertainty, beyond those from l/rrib expansion or QCD perturbation theory, 
have to be investigated. 

It had been realized that there are, in principle, long-distance contributions to B — > X s 7 and 
B — > X s e + e~ related to the cc intermediate state, which go beyond QCD perturbation theory 
and which are also not described by nonperturbative higher order contributions in the 1 /m^ 
expansion. They originate from the quark level process b — > sec — > s'y and have previously 
been estimated using model calculations based on cc resonance exchange. For B — > X s 7 one 
may have b — > s^fr — > sj, where the \1/ (or a higher resonance) converts into an on-shell photon 
(q 2 = 0). This requires extrapolating the \1/ couplings to far off-shell values (from q 2 = M\ to 
q 2 = 0). From estimates perfomed along these lines effects of typically ~ 10% at the amplitude 
level were found in [11], |12| . Although the accuracy of this approach is difficult to quantify, a 
sizable uncertainty from this source could not strictly be excluded. 

In the case of B — > X s e + e~ f\, the intermediate cc resonances can be on-shell for appropriate 
values q 2 of the dilepton invariant mass. In this region of q 2 , the resonance contributions form 
a very large background to the flavor-changing neutral current signal. Experimental cuts are 
therefore necessary to remove this part of the g 2 -spectrum in the analysis of B — > X s e + e~, 
which is then restricted to the range of q 2 below (and in principle also above) the resonance 
region. Far from the resonances the situation with respect to long-distance effects is analogous 
to the case of B — ► X s 7. The tails of the resonances can extend into this region, however, the 
Breit-Wigner form usually employed in modeling the resonance peaks cannot be expected to 
give a correct description of the tail region away from the peak. Model-dependent treatments 
of the long-distance effects related to the tails of cc resonances have been discussed in the 



literature [Q, [13|, M , but again cannot be considered as fully conclusive. 

Recently, it has been proposed to treat in a model-independent way, employing a heavy quark 

expansion in inverse powers of the charm-quark mass, the nonperturbative contributions in 



B — > Xg'y related to the intermediate cc state |T5], [16], |T7|, This interesting approach, 
which leads to 0(k 2 Q CD /m 2 ) corrections, has the advantage of being well defined and sys- 
tematic. It avoids problems of double-counting of certain contributions, which are present 
in model descriptions involving both charm-quark and hadronic charm degrees of freedom 

^For simplicity we write B — > X s e + e~ throughout the paper with the understanding that our discussion is 
equally applicable to B — > X s /i + /i~. 



1 



simultaneously. Furthermore, for cc states far off-shell, the quark picture appears to be more 
appropriate than a description in terms of hadronic resonances ty', ...). Irrespective of the 
ultimate numerical accuracy that can be expected from an expansion in Kqcd/^c-, we believe 
that this method can still yield a useful order of magnitude estimate of the effect, which is at 
least complementary, and probably superior, to alternative model-dependent calculations. 
The main purpose of the present article is to study nonperturbative long-distance effects in 
B — > X s e + e~ in a model-independent way, using the methods previously applied to B — > X^. 
We will calculate the 0(h 2 Q CD /m 2 ) corrections to the dilepton invariant mass spectrum and to 
the forward-backward charge asymmetry in B — > X s e + e~ . We also obtain the 0{K^ CD /m 2 c ) 
correction to the branching ratio for B — > X s 7 as a special case of this analysis. Contrary 
to model estimates (involving cc resonances) both the magnitude and the sign of the long- 
distance effect are calculable and we find AB(B — > X S ^)/B{B — > X s ^) = +0.03. Our result 
agrees in magnitude with previous calculations |L5], [T7], [TB| . 

We remark that in comparison to the case of B — > X s ^ and B — > X s e + e~, long-distance 
effects in B — > X s vv are strongly suppressed by factors of C(m^/M|). Obviously the reason 
is that neutrinos do not couple to photons but only to weak bosons, whose contributions are 
likewise completely negligible for the long-distance effects in B — ► X s ^ and B — > X s e + e~. 
From these observations it is clear that B —>■ X s vv is a particularly clean channel from a 
theoretical point of view. 

The paper is organized as follows. After this Introduction, in section 2, we briefly review 
nonperturbative effects in B — > X s ^ and B — > X s e + e~ . We also discuss the s6-photon-gluon 
vertex, which will be of central importance for our subsequent calculations. Using the results 
of section 2 in the limiting case of an on-shell photon, we derive the 0(AQ CD /ml) correction 
to the B — > X s 7 rate in section 3. The main results of our paper, the 0(k 2 Q CD /m 2 ) effects on 
the dilepton invariant mass spectrum and the forward-backward asymmetry in B — > X s e + e~ , 
are calculated and discussed in section 4. Section 5 contains some comments on long-distance 
effects in B — ► X s vv. A short summary is presented in section 6. 



2 Nonperturbative Effects in B —> X s 7* and the sb^fg— 
Vertex 

As mentioned before, B — > X s ^ and B — > X s e + e~ are dominated by perturbatively calculable 
short-distance contributions, including those from virtual top quarks. At NLO in QCD the 
intrinsic uncertainty of the perturbative calculation, indicated by residual renormalization 
scale dependence, is about ±6% for B — > X s j |T], ^ and of similar magnitude (or somewhat 
smaller) for the g 2 -spectrum in B —>■ X s e + e~ It is useful to keep these numbers in 
mind for comparison with uncertainties from other sources. Several effects that are beyond a 
perturbative treatment lead to corrections to the purely partonic decay picture and have to 
be considered for a more complete assessment of the total theoretical uncertainty: 

• Power corrections in l/m^: Sub leading corrections in the heavy quark expansion start 
at 0(k 2 Q CD /m 2 ) and amount typically to several percent. The effects are calculable, 
generally well under control, and can be taken into account in a complete analysis 0, 0. 



An exception is the photon energy spectrum in B — > X s ^ |L9| and the endpoint region of 
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Figure 1: The sb^g vertex at lowest order in QCD. 

the dilepton invariant mass (q 2 ) spectrum in B — > X s e + e~ J7|, where the 1/rrib expansion 
breaks down. This is not a problem for the integrated total B — > X s ^ rate and for the 
B — > X s e + e~ spectrum in the lower q 2 region. It is more troublesome for the shape of 
the photon energy spectrum in B — > X s 7, which is needed to extract BiB — > X s ^{) from 
the data. At present this fact still introduces some amount of model dependence in the 
measurement of this quantity [pOR . 

• Onshell cc resonances: In the case of B — > Xg'j there is a background from the process 
B — > X s ty, with the (on-shell) ^ subsequently decaying into a photon plus light hadrons. 
However, the energy spectrum of photons from this cascade process peaks at lower 
values than the one from the short-distance decay B — > X s 7. In the window 2.2 GeV 
< E y < 2.7 GeV used to measure B(B — ► X s ^) this background is fairly small 
and can be corrected for. In the case of B — > X s e + e~ the cc resonance background 
shows up as large peaks in the dilepton invariant mass spectrum. As discussed in the 
Introduction, it has to be removed by appropriate cuts. 

• Off-shell cc nonperturbative effects: A particular nonperturbative contribution to the 
short-distance process B — > X s 7 comes from the diagrams shown in Fig. [l|. It arises 
when the gluon is soft and interacts with the spectator cloud of the b quark inside the 
B meson. This is the effect pointed out in |y| |16[ and further discussed in |17], pL8]| . The 
same mechanism affects also B X s e + e~ , in which case the photon invariant mass is 
gV 0. 

It is this last point we would like to focus on in the present paper. Since the gluon is very 
soft and interacting with the light spectator cloud, the process in Fig. [I] corresponds to a cc 
pair from b decay converting into a (on- or off-shell) photon. Previous attempts to describe 
long-distance contributions based on the conversion of intermediate (cc) resonances into a 
photon (for both B — > X s 7 and B — > X s e + e~) can be viewed as model calculations aimed at 
addressing essentially the same effect. 

Staying with the quark picture, one may evaluate the diagrams in Fig. [I] directly. To first 
order in the gluon momentum k (with k 2 = 0), the vertex function can be written as 

A x sfng = -i^cC 2 egQ c sj,(l - l5 )G as b [e a ^%q x - e aX ^q 2 + e^% q a } . (1) 
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Here A is the Lorentz index of the photon, q the photon momentum and G aQ = G^ g T a the 
gluon field strength. C2 is the Wilson coefficient of the operator Q 2 = (sc)v-A(cb)v-A in the 
standard effective weak Hamiltonian. In the leading logarithmic approximation, 
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(2) 



Furthermore, A c = V* s V cb and Q c = +2/3 is the charm-quark electric charge. The form factor 
F(r) as a function of 

„2 

(3) 



Am 2 



is given by 



F(r) = — 



1 



arctan 



1 — r 



1 - Jl - l/r 
In V , + /- I - I 



< r < 1 , 



r > 1 . 



(4) 



2^r(r - 1) V 1 + V 1 - l / r 
Our sign conventions are such that e 0123 = +1 and the covariant derivative is 



D ll = d ll -igCr a Al + ieQ f A i 



(5) 



The diagrams in Fig. [I] generate a finite amplitude. However, in addition to the gauge- 
invariant m c -dependent structure in ([I]), also terms proportional to e aXtJ-e (q — k) a appear. 
We have not included these terms in A^ big because they are cancelled by the top-quark 
contribution via the GIM mechanism. This can be easily seen by considering the hierarchy 
m c <C m t Myy, which is sufficient since we are not interested in the detailed form of the 
top-quark contribution to A^ big and neglect 0(l/m 2 ) terms anyway. We note also that by 
virtue of the GIM cancellation between top- and charm-quark contributions, the integrals are 
all manifestly convergent and the full calculation, including the Dirac algebra, can be safely 
performed in four dimensions. 

Eq. (Jl]) will be sufficient for our actual calculations of 0(AQ CD /m 2 ) effects in B — * X s e + e~ 
and B —> X s ^. It is, however, useful to also consider the more general result for the stryg- 
vertex obtained without expanding in k ■ q. This will allow us to get some idea about the 
reliability of the lowest order approximation and its stability against higher order effects. Of 
course, there are further contributions, beyond those from higher powers in k ■ q. For instance, 
one could have terms with more than one gluon field. But these are suppressed by additional 
powers of AQ CD /ml and we shall neglect them. The vertex function (valid to all orders in q 2 
and k ■ q, but with k 2 = 0) is then given by 



A. 



sb^yg 



-^A c C 2 esQ c s 7M (l-75)T ffl & 1 
x F(r + t, t){G a Xa dPF ea + F^dPG* + G a d a F aX ) + F(r, -t)G a aQ d°F a 



(6) 
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where ^ 

2ml ^ 
and is the electromagnetic field strength. The form factor F is given byf] 
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The following relations hold 

F(r, 0) = F(r) , F(0) = 1 . (9) 

Expanding F in powers of t, can be understood as a sum over an infinite number of 



operators involving all powers of the gluon momentum [[L7|, |TJ| . Unfortunately, for the physical 
processes we want to study, only the lowest order contribution can be calculated. The higher 
derivative operators lead to matrix elements that are unknown and one has to resort to a 
merely qualitative discussion in this case. 

Since q is of order m;, and we are interested in the limit where k ~ O(Aqcd), we expect the 
matrix elements of these operators to be suppressed by corresponding powers of mbAQCD/fnl 
(suppression of t n+1 with respect to t n ) and Aq C D/ m b (suppression of FdG versus GdF). If 
both of the following inequalities were satisfied 

m b > Aqcd , (10) 
m 2 c /m b > Aqcd , (11) 

we could neglect all the operators involving higher powers of the gluon momentum and use (|Ij) 
instead of (f|). However, whereas ( |T0D is well justified, ( |TTD is not a priori a good approximation 
in the real world. Thus we can safely neglect the FdG term in (|6|) but we must have a closer 
look at the 0(t n ) corrections to (|T]). 

In the B — > X s ^ case (q 2 = 0) these corrections are estimated to be small [17], |18| . Indeed, 
even though k ■ q/m 2 c ~ 0(1), the Taylor expansion of F(t,t) in terms of k ■ q/m 2 c involves 
small coefficients 

rM=1+ *^ + 3(^y + «(^y + » (^y + ... (12) 

15 m l c 35 \ m z c J 525 \ m l c j 693 \ m z c J 



A typical value for t would be t m 0.3 [17|]. Then we have F(0.3,0.3) ~ 1.2, compared with 
F(0, 0) = 1. This indicates that the summed contribution of the t n -terms (n > 1) is probably 
not too important, owing to the small numerical coefficients in the series [|17|, Of course, 
this is a rather crude argument, since the matrix elements of the higher derivative operators 
may be weighted differently for different n. Furthermore, coherent addition of all the terms 
may not be correct either, but is at least a conservative assumption. (Indeed, the matrix 
element of the n = 1 operator turns out to be vanishing in the case of B — >• fT8||.) Still 



* Note that the form factor F can be written as F(r, t) — — 3847r 2 m^C , 2o('? 2 , k-q) in terms of the subtracted 
three-point function C20 defined in p2[ ; the explicit expression of C20 can be found in the second ref. in [E2[ . 
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this estimate, despite being rough, provides a certain measure of the possible importance of 
neglected effects, exploiting the information contained in the form factor. 

In the same spirit one may investigate the general case q 2 7^ 0. The ratio F(r+t, t)/F(r, 0), 
measuring the sensitivity to a nonvanishing t, reads ~ 1.3, 1.4, 1.5, 1.75 for r = 0.3, 0.4, 0.5 
and 0.6, respectively, using t = 0.3 (these values of r correspond to q 2 /ml ~ 0.1, 0.14, 0.17 and 
0.2). This suggests that the expansion in t makes sense provided that q 2 is far enough from the 
Am 2 threshold. On the other hand, due to the physical cuts at q 2 = Am 2 and {k + q) 2 = Am 2 
of the diagrams in Fig. [I], the functions F diverge near the 4m 2 threshold. This breakdown 
of the gluon-momentum expansion can be viewed as an indication of the appearance of large 
genuine long-distance effects due to the nearby cc resonances. Hence, we conclude that if 

Am 2 c -q 2 >m\ (13) 

the lowest order it = 0) result is still reasonable, as a sensible approximation for smaller q 2 
and as an order-of-magnitude estimate closer to this bound. 

So far we have been considering the region of q 2 below the resonance domain. Above the 
resonances r is parametrically large, r ~ q 2 /m 2 « m\jm 2 c 3> 1, and one may expand F(r+t, t) 
not only in t but subsequently also in inverse powers of r. Keeping only the leading term in 
1/r in each of the coefficients of t n one has 

OO O O J. 

Fir + t,t) « Y (-l) n+1 7 ^—rt n = -— + -! + ... (14) 

This shows that the leading corrections to the t = result behave as a series in powers of 
t/r ~ k ■ q/q 2 ~ KqcD/ m b- Numerically, using the full expression for F(r + t,t) and taking 
r = 2 as a relevant example, we find that F(2 — 0.3,-0.3) = —1.40 + l.lOi. To be more 
conservative we have here chosen t < so that all contributions add coherently. The result 
is reasonably close to the t = value F(2, 0) = —1.22 + 0.83i. Thus it appears justified to 
neglect the higher dimensional operators in this case and we expect the approximation t = 
to work in the large q 2 region of the spectrum. 

In the previous discussion we have restricted our attention to the function F(r +t,t), but 
similar conclusions hold also for F(r,—t). The latter contributes only for q 2 7^ and its 
expansion in t shows a better behavior than the one of F(r + t,t). 

At this point it is natural to address the question of the effect generated by the diagrams 
in Fig. [I] if the charm quark inside the loop is replaced by an up quark. At first sight, the 
explicit ml dependence of (^) would suggest that the effect diverges as l/m 2 when m q — > 0. 
This is still not a problem for B — > X^*, where the up-quark contribution is suppressed by 
the CKM hierarchy, but could be a dramatic effect in B — >• X^*. However, it should be 
noted that the l/m 2 behavior of (|) is correct only if the expansion of F in terms of k ■ q/m 2 
is allowed, which is not true for the up-quark contribution. In this case, and considering 
q 2 = for the moment, it is more appropriate to expand F(t,t) in inverse powers of t: 
the leading term in this expansion goes like t -1 ~ m 2 /k ■ q and cancels the artificial l/m 2 
divergence of (|6]). The operators generated by this expansion are nonlocal and their matrix 
elements cannot be estimated reliably. However, from naive dimensional counting, we expect 
the leading contribution to be of order Agc^/mf,. 

For the up-quark contribution in B — > X s ^e + e~ the situation is, in a sense, even more 
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favorable. If q 2 ^ (2GeV) 2 , one is above the region of uu resonances and q 2 itself provides 
the relevant short-distance scale. The leading long-distance contribution is then of the order 
K 2 Q CD /q 2 . Corrections to this behavior arise as powers of t/r ~ k ■ q/q 2 ~ Aq C d/ Vq 2 , 
corresponding to a series of local, higher-dimensional operators. These results follow from 
(14). The situation is analogous to the one for the high-g 2 end of the spectrum in the case of 
intermediate cc, except that for the uu case the operator product expansion approach remains 
valid down to smaller q 2 due to the absence of heavy resonances. 

We finally remark that the diagrams in Fig. [I] actually contribute in two distinct ways to 
the decays B — > X s j* . First, the gluon can be soft and couple to the light cloud in the B 
meson. This is the effect we are mainly interested in here and which we shall calculate, far 
from the resonance region, using ([]]). Second, the gluon, now not necessarily soft, may be 
radiated from the charm-quark loop to end up in the final state X s . This process is calculable 
in perturbation theory, it is infrared finite and contributes as part of the matrix element 
calculation at NLO in B — > [p], p3| . For B — > X s e + e~ the situation is analogous, except 
that the diagrams of Fig. |T| would enter only beyond the next-to-leading order (at 'next- 
to-next-to-leading order'). The first of these processes is not contained within the second 
one, which is entirely perturbative. Therefore in calculating the nonperturbative correction 
there is no double counting of contributions already taken into account in the higher order 
perturbative result. 



3 B X sl 

In this section we re-derive the 0(AQ CD /ml) correction to T(B — > X s ^) using the results of 
the previous chapter in the limit of an on-shell photon. 

By means of the optical theorem the total rate of a B meson decay can be expressed as the 
expectation value 

r=^-(B|T| B ) (15) 

of a transition operator T defined by 

T = Im i J d 4 x TH e ff{x)H ef f{0) . (16) 

Here 7i e // is the low energy effective Hamiltonian governing the decay under consideration. 
The states \B) in (|15D are to be taken in conventional relativistic normalization ((B\B) = 
2EV). 

Since we are analyzing a small correction to the full decay width, it is sufficient to work to 
leading logarithmic accuracy and to neglect consistently all relative 0(a s ) effects. In this 
approximation we may write 

Keff = ftfff + ft { eff ^ ( 17 ) 

where 

H ( 3 f = -%V: s V tb C 7 7 + h.c. , O r = ^-m b m^{l + l5 )b F, u , (18) 



represents the well known leading-log effective Hamiltonian for b —* 57 decay |24j (for a review 
see P5| ). C-j ~ —0.30 is the scheme independent ('effective') Wilson coefficient, whose analytic 
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form can be found for instance in 125 . 



The correction Tl^L can be derived from the amplitude (|I]) in the limit where the photon is 
put on-shell. In this case the vertex (|IJ) is equivalent to a local operator. Performing the 
usual matching procedure one obtains 



W$f = ^§Ok C 2 On + h.c. , (19) 



°n = 7TT-**y*0- ~ -Y^gGuxb e^d x F sa . (20) 



To leading order only Tie// contributes to (|15|), ([16|). One finds the well-known result (V* s V c b — 

-v,',v a ) 

T(B - x ^ = l5 (l/ ^ )2 V c? ' (21) 

The interference between Tr^jt and Ti-^ff leads to a correction of the transition operator 

AT = -w ( ^ )2 s%?^' G6 - = ■ < 22 > 



c 



Using 



we finally obtain 



(£|&0<7-G6|fi) 3 



2M B 2 



-(M^-Mi), (23) 



ffb| = -s^^r^ h +0 - 03 • < 24 > 

This result suggests that the potentially problematic nonperturbative effects in B — > X s 7 
should indeed be negligible, re-enforcing the role of the B — > X s ~/ process as a significant test 
of the Standard Model. 

We stress that the correction (23) has a definite sign, contrary to previous estimates of non- 



perturbative cc effects based on off-shell resonance exchange [|ll], [12| . In this context we note 
that for a positive matrix element (^3]) and negative CV, the consistent covariant derivative, 
specifying the sign of strong and electromagnetic couplings and thus of (f24|), is the one given 
in©. 

The nonperturbative correction in ( p^j ) may be compared with model calculations of the 
long-distance effect. Assuming that the dominant contribution to (p4|) is generated by the ^ 
exchange, the above result can be used to fix sign and magnitude of the nonfactorizable bs^fr 
coupling at q 2 = (g 2 in the notation of |12[). The result thus obtained is consistent with the 
lower values given in JT2] (g 2 ~ 10 



-2\ 



4 D — X s e+e" 

In the present section we proceed to compute the nonperturbative correction to the rare decay 
B — > X s e + e~. The effective Hamiltonian for b — > se + e~ at next-to-leading order has been 
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reviewed in [25]. In the following we will make use of the results collected in this article and 
adopt its notation. At next-to-leading order, the amplitude for B — > X s e + e~ can be written 

as 



A 



y/2 2vr 



2m b C 7 J™-^ wy u v + Cg ff Jg uj v v + C 10 Jg wy u j 5 v 



where 



<Jq 



(X s \sa w (l+ l5 )b\B) 
(X a \trf(l-y 5 )b\B) 



(25) 

(26) 
(27) 



and u, v denote the electron and positron spinor, respectively. C 7 and C w are Wilson coef- 
ficients; is not strictly speaking a Wilson coefficient, since it includes also contributions 
from the b — > se + e~ matrix elements of the operators Q±, ... Qe in the Hamiltonian p5 [. 
However, the amplitude expressed in terms of Cg provides a convenient notation. In partic- 
ular, Cg^ already contains contributions from intermediate charm-quark states entering the 
one-loop matrix elements (these are similar to the graphs in Fig. [l] but without the gluon 
and a lepton line connected to the photon propagator). 

In writing (p5|) we have used V t *V t b ~ — V* s V c b = — A c . Note also that, to calculate B — > X s e + e~ 
at NLO, only the leading logarithmic approximation is necessary for C 7 (= C^ eff in the no- 
tation of 

The nonperturbative correction we want to evaluate stems from the vertex function in (H), 
which yields the following correction to the leading amplitude in (E5j) 



\/2 



G F . a C 2 Q C F{r) 



'2ir 3m 2 . 



J, 



Q 



e ^ 8q 2 



wy\v 



(2* 



with 



Ja,^ e = (X s \s llM {l - lb )gG ae b\B) . (29) 

Note that the term proportional to q x in ([!]) vanishes when multiplied by uy\v due to current 
conservation. 

Adding A + A c , squaring the amplitude, summing over inclusive states X s and performing 
the necessary phase space integration, one may calculate the differential decay rate for B — »■ 
X s e + e~. Defining s = q 2 /m 2 and 



R(s) 



d T(B 



ds 



X q e + e- 



T(B -> X c ev) 

one has at NLO in QCD perturbation theory [z = m c jm^) 

fl W^7(^S[ (1 + 2s) ( |<5 » //|2 + l<5 '° |2 ) +4 ( 1 + 

Here f(z) is the phase space factor and k(z) the QCD correction factor entering T(B 



\C 7 



2 + 12C 7 ReCg ff 



(30) 

• (31) 
X r ev) ; 



they can be found in [25 



For the 0(h. 2 Q CD /m 2 ) correction term, arising from the interference of A c with A, we obtain 



AR(s) 



or „ 2(Af|, - M 2 B ) (1 - sY 



4tt 2 



9m;? 



Re \ F(r) 



/(*)"(*) I 



a 



1 + 6s - s z 



+ C 9 e// *(2 + s) 



(32) 
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Figure 2: The dilepton invariant mass spectrum of B — > X s e + e~ normalized to the semilep- 
tonic width ([!(]). The solid curve is the short-distance result ([!!]) whereas the dashed one 
includes the 0(k^ CD /m 2 c ) corrections. The dotted lines indicate the position of the \1/ and \P' 
resonances. A dash-dotted line is shown at q 2 = 3m 2 . The corrections diverge at q 2 = 4m 2 . 
In this plot we have used m c = 1.4 GeV and m b = 4.8 GeV. 



For the evaluation of AR we have used the identity |26| 

(B\lYG aP b\B) = — (B\ba ■ Gb\B) tr {r(l + 
48 



cr. 



a/3 



(33) 



which is valid to leading order in the 1/m;, expansion. Here v denotes the B meson four- 
velocity and T an arbitrary string of 7-matrices. For handling the Dirac algebra throughout 



our calculations we have used the program Tracer [[27]. 

The final result is displayed in Fig. ^ showing the corrected and uncorrected R(s). The 
conditions for the validity of the corrections have already been discussed in section 2. They 
are best satisfied for the high-g 2 end of the spectrum, above the resonances, and for very low 
q 2 . The reliability of the approximation used deteriorates somewhat as q 2 is increased from 
~ towards the resonance region. As an order of magnitude estimate the calculation should 
still make sense around q 2 w 3m 2 , indicated in Fig. ^ by the dash-dotted line. For larger 
q 2 one gets too close to the lowest resonance and the heavy quark expansion breaks down. 
Further up, at q 2 = 4m 2 , the correction develops a (unphysical) square root divergence (in our 
choice m c = 1.4 GeV). The difference between and 2m c is of 0(Aq C d) and corresponds 
to nonperturbative effects that are beyond the control of our approximation. 
As can be noticed, the 0(h 2 Q CD /m 2 ) corrections are very small in the region where this calcu- 
lation can be trusted, that is below 3m 2 and above the resonance peaks. Actually, in the high 
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q 2 region the effect is further suppressed since F(r) decreases as 1/r for large r, as discussed 
in section 2. For q 2 ~ 0(m 2 ) the corrections behave in fact more as 0(h 2 Q CD /m 2 ) instead of 
0(AQ CD /ml). We also observe that for s below 0.1 the corrections become vanishingly small 
due to a cancellation among the terms in (|3^) (CV is negative, KeCg^ positive). This is related 
to the fact that the correction AR(s) is negative for s above 0.1, where Cg dominates, but 
positive for s close to zero, which essentially corresponds to the case of B — > X^. In fact, 
one may obtain (g3|) special case of AR/ R from fl3~l|) and (^) in the limit q 2 — > 0. 



The sign of the effect for s around 0.2 is different from the one obtained by the approaches 
based on resonance exchange 0, 13, 14]. This is not a problem since the sign of the correc- 



tion determined in model calculations is not reliable far from the resonances. Moreover, as 
discussed in the previous section, the corrections we are considering are essentially related to 
the nonfactorizable contributions of charmed resonances. The latter are usually considered in 
B — > X s j, where the factorizable terms vanish, but are always neglected in B — > X s e + e~. 
The factorizable resonance contributions can be identified, to the lowest order in a s , with 
diagrams as those in Fig. [I] but without gluons. The effect of such diagrams is positive and 
is already included in the NLO calculation as the matrix element contribution to Cg . As 
pointed out in P|| within the context of exclusive decays, adding the factorizable resonance 
effects to the NLO calculation leads to a double counting problem related to the simultaneous 
use of quark and hadronic (cc resonances) degrees of freedom. On the other hand, this is 
not the case for the Ag CD /m^ effect we have evaluated, where we have consistently used a 
pure quark level description. Of course, some model-dependent treatment may still be useful 



close to the \l/-resonance (an interesting approach has been presented in [^9]). If one should 
attempt this, the 0(k 2 Q CD /m 2 ) effect could be used to fix the small q 2 limit. 

Another interesting quantity that can be measured in B — > X s e + e~ decays is the forward- 
backward charge asymmetry [ 13| . Normalizing T(B — > X s e + e~) to the semileptonic width, as 



in (|30|), we define 



A{s) = — — — / dcosO 1 s ^sgn(cos^), (34) 

1 [B — > X c ev) J -i as a cost) 

where 6 is the angle between the e + and B momenta in the dilepton center-of-mass frame. 
Proceeding similarly to the case of R(s), the NLO perturbative result turns out to be 

whereas the 0(AQ CD /ml) correction term, arising from the interference of A c and A, is given 
by 

A ^-g c - j %^7^ Re ^ F(r) } (1 + 3s) ' (36) 

The corrected and uncorrected results for A(s) are shown in Fig. [| Concerning range of 
validity and size of the corrections similar comments apply to A(s) as to the spectrum R(s) 
discussed above. In addition we note that the small corrections to A(s) and R(s) have a 
tendency to further cancel in the normalized asymmetry A(s)/R(s), a useful observable for 
direct comparison with experiment. 
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Figure 3: The differential forward-backward asymmetry defined in (|3lf) as a function of the 
dilepton invariant mass. Notations are as in Fig. |2|. 

5 B -> X s vv 

In B — ► X^z/ there are no contributions from virtual photons. Consequently the pattern of 
GIM cancellation is not logarithmic ('soft'), as in B — > X s e + e~, but powerlike ('hard') due 
to the exchange of heavy gauge bosons. Therefore it is clear that long-distance effects in 
B — ► X s uu from charm quarks are additionally suppressed by a factor of order vr? c jMy V . In 
principle, however, also in this case long-distance contributions related to the cc intermediate 
state do exist. The largest of these is associated with the cascade process B — > X s ty — > X s vv. 
In the following we shall try to estimate this effect to show more quantitatively to which 
extent the short-distance contribution dominates in B — > X s vv. 

Using the relation (0|c7 M c|\l/) = m^f^e^ and the coupling of the cry^c current with the elec- 
troweak fields 



£ = —ecj^c 



2 * 1 

3 M 4 sin 9,„ cos 9„ 



sin 9,i 



we obtain 



vv) 



fM 9 



3-8 sin 2 9„ 



\A(*^e+e-)\ \M Z 

Then, using the theoretical result B s d (B 
mental values of B(B -> X s ^) and B(V -> e+e~) 0, we find 

,B(B -> X s ^)5(# 



3.3 x 10" 



(37) 



(3? 



16^2 cos 2 6L sin 0„ 

X s z/z/) = (4.0 ± 1.0) x 10" 5 § and the experi- 



B(B -> X^(z/z/)) 
B s d (B -> X s i/z7) 



3-R,, 



e 1 e 



B sd (B 



X s vv) 



5 x 10" 



(39) 
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Unfortunately we cannot reliably estimate the interference terms between short- and long- 
distance contributions, since we have no useful information about the strong phases in B — > 
X s ty. However, the result ( p9| ) indicates that the long-distance corrections to the amplitude 
are at the level of 10~ 3 at most. This is by far enough to conclude that B — > X s vv is an 
extraordinarily clean channel || and, therefore, a very interesting probe of the Standard Model 
and its extensions |§. 

6 Summary 

In this paper we have calculated long-distance corrections of 0(k^ CD /m 2 c ) to the dilepton 
invariant mass (q 2 ) spectrum and to the forward-backward asymmetry in B — > X s e + e~ . For 
low q 2 jm\ ^ 0.1 these corrections can be expected to give a fairly accurate description of 
long-distance effects due to intermediate (far off-shell) cc pairs. Contributions from higher- 
dimensional operators, formally suppressed by powers of m^Ag^^/m;; relative to the XQ CD /m 2 
effect, involve unknown hadronic matrix elements and have been neglected. Despite the fact 
that mbAQCD/m 2 ~ 0.6, this is justified for small q 2 due to small numerical coefficients accom- 
panying these additional contributions. They could be more important for larger q 2 . Close 
to the resonance region q 2 /m 2 ~ Am 2 /ml ~ 0.34 the assumptions on which the calculation is 
based clearly break down. On the other hand, simple order-of-magnitude estimates suggest 
that even for q 2 /m\ ~ 0.26 the calculable leading 0{Is?Q CD /m 2 ) effect should give a reason- 
able account of the approximate size of the corrections within, say, a factor of two. Above the 
resonance region q 2 /m 2 ^ 0.7 the calculated effect is again reliable. In spite of the various un- 
certainties we believe that the l/m c expansion approach is still preferable to model-dependent 
estimates relying on off-shell cc resonances. 

Numerically the corrections we find are small (at the one to two percent level) in the 
phenomenologically interesting region of q 2 away from the resonances. For comparison one 
may recall that uncertainties from QCD perturbation theory are at the level of ±6% at next- 
to-leading order. This indicates that long-distance corrections related to intermediate cc 
should not be a serious problem in the analysis of B — > X s e + e~ decay, unless the neglected 
effects of higher order in AQ CD /m 2 or m^QCDj'^'c turned out to be substantially larger than 
anticipated. 

We have pointed out that the 0(k 2 Q CD /m 2 ) effect is calculable in both magnitude and sign, 
which is to be contrasted with the situation in hadronic model calculations. The 0(XQ CD /m 2 ) 
correction to B(B — > X s j), considered in previous work, can be obtained as a special appli- 
cation of our analysis and is found to be ~ +3%. 

The decays B — > X s 7 and B — > X s e + e~ may be contrasted with the mode B — > X s vv, 
where the corresponding nonperturbative effects due to intermediate cc are still much stronger 
suppressed, down to a level of ~ 10 -3 in the amplitude. 

Long-distance effects are in general a notoriously difficult problem. However, they have to 
be controlled sufficiently well to address fundamental questions in flavor physics. The existence 
of a realistic limit in which such effects can in fact be computed from first principles is therefore 
of considerable theoretical interest in its own right, as well as of interest for phenomenology. 
It is particularly gratifying that for the important processes B — > A s 7 and B — > X s e + e~ a 
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calculation of this type is possible and a class of potentially large long-distance corrections is 
indeed found to be rather well under control. 

Note added: After this paper was finished, we have received a preprint |Hj] in which the 



0{K^ CD /m 2 c ) correction to the differential decay rate of B — > X s e + e has been calculated. 
We disagree with their result in sign. 
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